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The question of thermodynamic consistence and $-derivability of the asymptotic limit of high 
spatial dimensions for quantum itinerant models is addressed. It is shown that although the irre- 
ducible n-particle Green functions are local, reducible vertex functions retain different momentum 
dependence. As a consequence, the vertex corrections to conductivity do not generally vanish in the 
mean-field limit. The mean-field theory is a <I>-derivable approximation only if regular nonlocal or 
anomalous local external sources are admitted. 
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The Baym and Kadanoff concept of thermodynamic 
consistence and conserving character of approximations 
is derived from an observation that knowing a thermo- 
dynamic potential as a functional of external distur- 
bances we can reconstruct thermally averaged values of 
any observable [Q. The most transparent way to guar- 
antee thermodynamic consistence is to find a generating 
Luttinger-Ward functional $ [Gj] depending on the fully 
renormalized one-particle propagator Gj in the presence 
of a set of external disturbances J a . Functional deriva- 
tives of the generating functional with respect to G lead 
to irreducible functions, while functional derivatives with 
respect to the external sources J a to reducible ones 0. 
The n-th orders of the derivatives generate n-particle 
functions, either irreducible or reducible. This definition 
of higher-order Green functions guarantees a conserving 
character of the approximation. But it may be in conflict 
with criteria for selecting Feynman diagrams contribut- 
ing to the n-particle Green functions in special limits 
such as the mean-field theory, ft is defined as an exact 
solution in infinite spatial dimensions. 

Selection of relevant diagrams in high spatial dimen- 
sions is determined by a scaling of nonlocal terms in lat- 
tice Hamiltonians. It can, from principal reasons, be 
made only within a perturbation expansion. In lattice 
models of correlated and/or disordered electrons one can 
conclude from finiteness of the energy density on a col- 
lapse of diagrams to a single site, || or on irrelevance 
of momentum conservation at internal vertices |Q. The 
one-electron self-energy gets local in infinite spatial di- 
mensions and can be generated from a Luttinger-Ward 
functional depending on only local elements of the full 
one-particle propagator. Such a functional, or its Legen- 
dre transform was constructed for Hubbard-like models 
without and with chemical randomness in H[|. As far 
as we are interested only in the dynamics of one-particle 
quantities, we can completely neglect the off-diagonal el- 
ements of the propagator |Q . 

The limit of high spatial dimensions has an asymp- 
totic character, i. e. leading orders in the small parame- 
ter dT 1 are to be kept for each quantity. The off-diagonal 
elements of the one-particle propagator scale on a hyper- 



cubic lattice as gM^I/ 2 , where |i— j| the lattice distance 
between sites i and j. The off-diagonal elements are es- 
sential for susceptibilities deciding about the existence of 
solutions with a long-range order ||. 

To prove a long-range order we have to evaluate two- 
particle susceptibilities (correlation functions) at fixed 
momenta. The collapse of diagrams in the mean-field 
limit is well defined in the lattice space but the cor- 
responding reduction in the momentum dependence of 
two-particle Green functions, depending generally on 
three momenta, is less evident. Particularly, if the limit 
d — > oo is used within the Baym and Kadanoff conserv- 
ing scheme, the locality in two-particle functions is still 
understood controversially pJlQ||. 

It is the aim of this Letter to address the problem of 
compatibility of the high-dimensional asymptotics with 
the Baym concept of $-derivability of all Green functions. 
We show that all nonperturbative irreducible vertex func- 
tions are local and can be generated from the Luttinger- 
Ward functional via derivatives w.r.t. local propagators. 
Reducible n-particle Green functions are generally non- 
local, but the nonlocality is perturbative. They can be 
derived from the generating functional if we either keep 
the external sources nonlocal or allow for local anomalous 
perturbations that do not conserve charge or spin. 

We can construct the limit to infinite lattice dimen- 
sions within the Baym and Kadanoff scheme. We start 
with a general model the partition function of which can 
be represented via a functional integral 

Z{J;G(°)- 1 }= (1) 

VtpVtp* exp {-93*nG (0) - V + <P*J<P + U [<p* , <p]\ 

where, for the sake of simplicity, we suppressed all deco- 
ration indices for the lattice and internal degrees of free- 
dom. All the interaction is contained in a local functional 
U , J stands for external sources, and ?y = ±1 for bosons, 
fcrmions, respectively. 

A renormalized perturbation theory is achieved if 
we replace the bare one-particle propagator G^ with 
Q = (G^ 1 + It is a straightforward task to write 
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down a generating functional for the renormalized pertur- 
bation theory. The grand potential contains the renor- 
malized one-particle propagator G and the self-energy E 
as variational functions as follows 

fi[J;G,E] = -/rMnZfOjG -1 + E} + V /3~hrlnG 

+r ? /r 1 trln [g^ 1 - 77 J - si . (2) 

The physical values of G and E are obtained from the 
equations <5<3>/ <5S — and <5$/(5G = 0, respectively. Note 
that G and S are Legendre conjugate variables. 

Representation (|l|) holds in any dimension. In the 
limit d — > oo we have to separate site-diagonal and off- 
diagonal parts of G and E. The electron systems show 
the following asymptotics in high dimensions 

G = G diag [d°] + G off [d-^ 2 } , E = E*° fl [d°] + E o// [d" 3 / 2 ] 

(3) 

The grand potential in the limit d = oo is obtained if 
the off-diagonal parts G°^ and E°^ are completely ne- 
glected in (|J) . The nonlocal part of the grand potential is 
contained solely in the bare propagator G' ' . Hence, all 
nonlocal quantities in the mean-field theory are treated 
perturbatively without renormalizations. 

According to Baym and Kadanoff, n-particle Green 
functions can be generated via successive applications of 
functional derivatives w.r.t. the external source J. Local 
sources are dominant in d — oo and we have for one- and 
two-particle functions 



G lv [J}=-(3 S -^- = (^ v )j, L n ',22< 
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(4) 



where 1 = (Ri,ti, ax, . . .) denotes a set of degrees of 
freedom in the space-time representation. The bar indi- 
cates a coordinate of a complex-conjugate field and any 
decoration of the labeling index refers always to the un- 
decorated lattice site. 

We now find equations of motion for two-particle func- 
tions and an explicit representation for the two-particle 
vertex in d = 00. As a first step we subtract the free 
two-particle propagation from L\v 22 1 an< i define a new 
vertex function 
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where the inversions use only the site-diagonal functions. 

Analogously we define a local vertex Jii' I"V" from 
an averaged cumulant two-particle function 



7ii\i"i'" = V 

2, 2', 2", 2 



We use the above definitions in the mean-field limit of the 
generating functional (||) with inhomogeneous but local 
one-particle functions to find an explicit form of a two- 
particle function bY,\v I needed for the evaluation 
of Lxx> 22'- If we switch from L to T we finally obtain a 
Bcthc-Salpctcr integral equation of motion 

E 1 8a,ifo,i'-v E 7iim<<i»<g^4g^, 1 r3 3 ,, 22 , 

3,3' [ J 

= <5ri,R 2 7I1',22' • (7) 

This form of the Bethe-Salpeter equation is natural for 
the mean- field limit, since it separates diagonal and off- 
diagonal propagators. The former are contained nonper- 
turbatively in the local vertex 7 while the latter partic- 
ipate only via two-particle bubbles forming the nonlocal 
part of the vertex T. 

We can further introduce a two-particle irreducible 
vertex that we denote A. This function enables us to 
replace in (^) the off-diagonal one-particle propagators 
with the unrestricted ones. A new equation of motion 
with the two-particle irreducible vertex reads 

^ \ <53,lfe',l' - I) E ^■lV,V'V"Gx'"iGyx» > rg 3 / ! 52' 
3,3' [ J 

= <5ri,R 2 Ai1',22' • (8) 

This form of the Bethe-Salpeter equation is common in 
general treatments without the collapse of diagrams. In 
the mean-field limit, d — 00, the irreducible vertex A is 
the same for 7 as well as for V . 

The irreducible vertex function is related to the self- 
energy via a generalized Ward identity 
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X <5r!,R 2 Gj^G^, {(^2>2'V4'V?2'''} 

{<f%<P2'){<pi»<P2'") - v{vi<P2"'){vi" ( P2')}Gj / } 2 , / G^j, 1 , /l . (6) 



following directly from the above derivation. Since the 
irreducible vertex functions appear in integral kernels of 
the Bethe-Salpeter equations of motion for higher-order 
Green functions, they get local in the lattice space in 
the mean-field limit and hence independent of momenta 
in spatially homogeneous solutions. Variations at fixed 
momenta as used in [p 10 loose their meaning where off- 
diagonal elements in the lattice space are treated pertur- 
batively. Nonlocal corrections to the irreducible vertex 
functions are negligible in leading order of d — ► 00 . 

The above derivation of the two-particle vertex func- 
tion within the Baym and Kadanoff conserving scheme is 
a generalization of the Brandt and Mielsch construction 
of charge susceptibilities in the Falicov-Kimball model 
S. The review on the d = 00 approach |Q uses, how- 
ever, a different construction of two-particle functions. 
The question is whether the two constructions lead to 
the same vertex function. In other words, does ^xi' 22' 



2 



contain all leading d^ 1 contributions to the two-particle 
vertex? 

It is undisputable that maximally two different lattice 
sites are relevant for two-particle functions in d = oo [|TT| . 
Two-particle Green functions have four end points. We 
generally have three possibilities how to select couples of 
end points from the same lattice site, Fig. [I]. 



J J 1 J J i 

FIG. 1. Leading order of two-particle vertex functions in 
d — » oo. The choices of two distant lattice sites i, j correspond 
to different two-particle reducibility channels, interaction (U), 
electron- hole (eh), and electron-electron (ee), respectively. 

Attaching end points of two-particle functions to lat- 
tice sites determines multiple repetition of nonlocal two- 
particle bubbles in the Bcthc-Salpctcr equations for two- 
particle reducible functions. Each choice in Fig. [I] is of 
the same order in the parameter d^ 1 and represents a 
different two-particle reducibility. The first choice, re- 
sulting from the above Baym and Kadanoff construction, 
leads to the so-called interaction two-particle channel of 
a chain of spin-triplet polarization bubbles screening the 
bare interaction. The second type of nonlocality in the 
two-particle function induces electron-hole, and the last 
one then electron-electron (spin-singlet) multiple scatter- 
ings. 

We use the three ways of constructing nonlocal two- 
particle functions and define three types of multiplication 
of the off-diagonal bubbles and the local vertex in the 
Bethe-Salpeter equation (0) 
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The three functions k" form three integral kernels of 
the Bethe-Salpeter equations for three vertex functions 
r Q . We hence have three different nonlocal vertex func- 
tions being of the same order in the small parameter d . 
These functions differ in their nonlocal parts and reduce 
to the same vertex 7 when all four end points are from 
the same site. The full vertex function is then a sum of 
the solutions to the three Bethe-Salpeter equations from 
which the local vertex 7 must be twice subtracted. 

It is more instructive to give an explicit representation 
for the full two-particle vertex function in the Fourier 
representation with fermionic k — (k, ito n ) and bosonic 



q = (q, iv m ) four-momenta. The three contributions to 
the full vertex function in A-space are obtained by Fourier 
transforming (^). Two-particle bubbles 

Xa a '{<¥,iu n ,iu nl ) = -^-^2G a (k,iuj n )G ljl (c[±k,iu} n ,) 



(11) 



carry the surviving momentum dependence of the two- 
particle functions in the infinite-dimensional limit. We 
have an explicit representation j|] for a hypercubic lattice 
with the nearest-neighbor hopping t = l/Vd 

/oo 
d\d\'8(\uj n )9(\'uj n >) 
-OO 

exp ji (Xx a (iuj n ) + A'av(iav)) - j (A 2 + A' 2 + 2AA'X(q))| 

(12) 

with x a (z) — z + /j, + aB — Yi a (z) and X(q) — 
5Zj/=i cos(q„). The last quantity, measuring nonlo- 
cality in infinite dimensions, behaves as a Gaussian ran- 
dom variable with variance d~ l when summed over mo- 
menta. The local bubble, entering the local vertex 7, is 
a product of local propagators and equals x with X = 0. 

ak "k + q 



a'k' 



a'k' + q 



FIG. 2. Generic vertex function in fc-space with three 
independent four-momenta and a defined order of incoming 
and outgoing fermions. 

To clarify the meaning of different four-momenta in 
the vertex functions we introduce a generic notation 
shown in Fig. ^. The full vertex function in leading 
asymptotic d~ 1 order in the absence of external sources 
J a reads in this notation 

r CTCT /(k, iw n ;k', iw n >; q, iv m ) = Y v aiJ , (q; iuj n ,iiv n >, iu m ) 
+T e * a , (k + k' + q; iuj n , iuj n ,,iv m ) - j a<T > (ioj n , iu> n > , iv m ) 
+r^,(k- ]s!;iu} n ,i(jj n >,iv m ) - ^ aal (iu} n ,iuj n >,ii/ m ) (13) 

where the each momentum-dependent contribution is de- 
termined from a separate Bethe-Salpeter equation (j^) 
with the appropriate multiplication from ([l0]). Each of 
the nonlocal terms in ( |l3| ) depends only on one transfer 
momentum reflecting the fact that maximally two differ- 
ent lattice sites are relevant. However, because of differ- 
ent ways to connect the end points in the two-particle 
vertex, the full function depends on all three indepen- 
dent momenta. We have 7(. . .) = TV -1 ^ q r Q (q; . . .) for 
all channels a = U,eh,ee. 
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To construct a Luttinger-Ward functional $ generat- 
ing higher-order functions of the mean-field theory, we 
have to account properly for the off-diagonal propaga- 
tors. Unlike the diagonal elements, suffering no restric- 
tion in d = oo, the off-diagonal propagators contribute 
only under favorite circumstances, i. e. when properly 
combined with summations over neighboring sites. That 
is why the off-diagonal elements appear in a perturbative 
manner in the mean-field limit. However, going to higher- 
order Green functions we have to take into account higher 
powers of the off-diagonal elements. The two-particle 
functions are generally of order d^ 1 . They can be gener- 
ated via functional derivatives of the local, but inhomo- 
geneous self-energy, or by derivatives of the nonlocal self- 
energy w.r.t. nonlocal sources or the off-diagonal propa- 
gator. The leading order of the nonlocal part of the self- 
energy is <i~ 3 / 2 and of the off-diagonal propagator d" 1 / 2 . 
Hence, a functional derivative JSG"** oc dr 1 con- 

tributes in leading order to the two-particle vertex. To 
turn the mean-field theory a ^-derivable approximation, 
we have to consider nonlocal external sources con- 
necting fluctuating fields at different lattice sites. The 
two-particle functions not derivable with local external 
sources can then be generated via nonlocal derivatives 
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Generally, n-particle functions demand to keep nonlocal 
sources in $ up to the nth power. Once we introduce 
nonlocal external sources into <i>, we have to expand all 
quantities, including the one-particle propagator and the 
self-energy appearing in $ in powers in J " . 

This perturbative construction with off-diagonal ex- 
ternal sources is applicable only in the high-temperature 
phase. At low temperatures, each constituent of the two- 
particle vertex function (13) controls a different type of 
fluctuations that may cause phase transitions. A singu- 
larity in the function T u signals a longitudinal, in T eh 
a transversal (magnetic) order, and a singularity in T ee 
indicates a superconducting instability. Neglecting any 
of the components of the full vertex function means sup- 
pressing potentially relevant fluctuations. 

The latter two transitions give rise to anomalous 
Green functions in the ordered phase. To describe these 
phases correctly we have to introduce anomalous (com- 
plex) local sources into the partition sum connecting, in 
a hermitian manner, </?|<^ for the transversal order and 
in the superconducting phase ]T^] . These anoma- 
lous perturbations induce local anomalous propagators 
and self-energies as variational functions in the gener- 
ating functional The mean- field Luttinger-Ward 
functional, that is able to generate the complete lead- 
ing asymptotics of higher-order Green functions and to 
describe all low-temperature phases, must contain either 
regular (of density type) nonlocal, or anomalous (com- 



plex) but local external sources and their conjugate vari- 
ational functions. However, only local variational func- 
tions and order parameters are treated nonperturbatively 
in d = oo. 

The newly derived vertex function ( filf ) contains the 
full momentum dependence and fulfills all the symmetry 
transformations of the exact solution. It means that the 
current-current correlation function and consequently the 
electrical conductivity do have nontrivial vertex correc- 
tions contrary to general expectations 0,0- It is only 
the function T eh that looses vertex corrections to the con- 
ductivity due to the symmetry. The other two, explic- 
itly depending on q, correct the conductivity of a single 
particle-hole bubble in high dimensions. 

To conclude, we showed how to construct a generat- 
ing functional for the mean-field theory with the correct 
asymptotics of higher-order Green and correlation func- 
tions. The mean-field theory as a complete solution in 
d = oo is a ^-derivable theory if either regular nonlocal 
or anomalous local sources are introduced. Only then we 
are sure to generate the complete leading asymptotics of 
higher-order vertex functions and to describe all possible 
low-temperature phases. However, the full higher-order 
vertex functions in d — oo cannot be generated from sin- 
gle equations or directly via functional derivatives from 
the Luttinger-Ward functional as in the exact theory. 
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